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ABSTRACT 


Near- optimal three dimensional trajectories from a low earth park orbit 

inclined at 28.5 deg to a synchronous -equatorial mission orbit have been 

developed for both the storable (thrust = 28,912 N (6,500 lbs), I = 

sp 

339 sec) and cryogenic (thrust = 44,480 N (10,000 lbs), I = 470 sec) 
Space Tug using the iterative cost function minimization technique con- 
tained within the Modularized Vehicle Simulation (MVS) Program. The 
finite burn tiroes, due to low thrust- to-weight ratios, and the associated 
gravity losses are accounted for in the trajectory simulation and optimi- 
zation. The use of an ascent phasing orbit to achieve burnout in synch- 
ronous orbit at any longitude is investigated. The ascent phasing orbit 
is found to offer the additional advantage of significantly reducing the 
overall delta velocity by splitting the low altitude bum into two parts 
and thereby reducing gravity losses. 
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INTRODUCTION 


One intended Space Tug mission which receives much attention due to the expected 

j 

frequency of use is the geosynchronous mission orbit. For this mission the 
Space Tug and attached payload are injected into a 28.5 deg inclined, low earth 
orbit by the Space Shuttle. The Space Tug then performs a maneuver to enter 
a transfer orbit with apogee at synchronous altitude. Upon reaching apogee, 

-K* 

the Tug's engine fires again to circularize the orbit and the payload is 
subsequently deployed. After a series of phasing maneuvers the Tug retrieves 
a second payload from synchronous orbit and performs a retrothrust deorbit 
burn to enter a transfer orbit with a low altitude perigee. At perigee the 
Tug burns to enter a phasing orbit which will produce the correct phasing 
relationship between the Tug and the waiting Shuttle for rendezvous purposes. 
Having completed a revolution in the phasing orbit, the Tug enters an orbit 
coelliptic with, and about l8.5 km (10 n mi) above that of the Shuttle, and 
acts as the passive vehicle in the ensuing coelliptic rendezvous. The geo- 
synchronous mission profile for the cryogenic Space Tug is shown in Fig. 1. 

The purpose of the present study is to determine the AV required by the Space 
Tug to perform the ascent portion of the geosynchronous mission using the 
minimum fuel trajectory. Trajectories are developed for both the cryogenic 
Space Tug and the storable Space Tug. Baselines for both these vehicles are 
contained in Refs. 1 and 2, respectively, and baseline data pertinent to this 
study is shown in Table 1. In both cases a shuttle weight constraint of 
29,484 kg (65,000 lbs) was assumed. This was also assumed to be the initial 
ignition weight for the Tug resulting in conservative losses. 

From Table 1 we can see that for both Space Tug configurations the thrust - to~ 
weight ratios are in the range where gravity losses due to finite burning times 
will be significant. For this reason the Modularized Vehicle Simulation (MVS) 
Program (Ref. 3) was used to accurately integrate the trajectories through 


This maneuver includes plane changes to achieve the geosynchronous- equatorial 
mission orbit. 














powered and coasting flight, thus including these gravity losses in computing 
the required AV. 

Table 1. Tug Baselines 

* 

Storable Tug Baseline 

Thrust - 28,912 N (6,500 lbs) 

Specific Impulse - 338 sec 

Shuttle Park Orbit - 278 x 278 km (150 x 150 n mi) 

i = 28.5 deg 

Mission Orbit - 35,787 x 35,787 km (19,323 x 19,323 n mi) 

i = 0 deg 

Cryogenic Tug Baseline 

Thrust - 44,480 W (10,000 lbs) 

Specific Impulse - 470 sec 

Shuttle Park Orbit - 296 x 296 km (l60 x 160 n mi) 

i = 28.5 deg 

Mission Orbit - 35,863 x 35,863 km (19,364 x 19,364 n mi) 

i = 0 deg 

* 

Taken from Ref. 2 

** 

Taken from Ref. 1 

In order to reduce these gravity losses, a three bum ascent, in addition to 
the nominal two bum ascent shown in Fig. 1, is analyzed. The three burn ascent 
which reduces gravity losses by splitting the large low altitude bum into two 
smaller burns which are performed closer to perigee, is shown schematically in 
Fig. 2. The first bum produces an intermediate orbit with a specified apogee 
altitude and accomplishes a small amount of plane change (generally of the 
order of one degree). Due to the finite burn time perigee altitude is also 
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Produces intermediate orbit with specified apogee altitude, incidentally 
increasing the perigee altitude. A small portion of the total plane change 
is accomplished . 

A burn straddling perigee produces a transfer orbit with apogee at synchro- 
nous altitude. Total plane change is increased to 2 deg. 

A burn initiated near apogee circularizes the orbit and accomplishes 
the remaining 26.5 deg of plane change. 


Figure 2. 


Three Burn Ascent to Synchronous -Equatorial Orbit 



raised. The second burn straddles perigee and raises apogee to synchronous 
altitude. The total plane change is increased to 2 deg. Slightly before 
reaching apogee, the circularization burn is performed and the remaining 
26.5 deg of plane change is accomplished. In addition to providing a lower 
total AV, the three bum mission has the additional advantage of offering an 
ascent phasing orbit which can be selected to achieve longitude phasing in 
synchronous orbit. Thus, optimal low thrust two and three burn ascent tra- 
jectories will be developed for both cryogenic and storable Space Tugs using 
the Modularized Vehicle Simulation Program. For the three burn ascent 
missions, the optimal intermediate orbit apogee altitude will be sought as 
part of this study. 
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2 . 


METHOD OF ANALYSIS 


The Modularized Vehicle Simulation Program has available an option to minimize 
a cost function specified by the user by means of an iterative convergence 
subroutine « For the purpose of this study the following cost function was 
specified for minimization: 


J = (A y ) + t«2 ( h a " h s ) + <*>3 ( e ) + (*) 


where 


* 


AV 

h 

a 

h 


s 

e 


i 


“J 


delta velocity expended in reaching final orbit 
apogee altitude of final orbit 

synchronous altitude 

eccentricity of final orbit 
inclination of final orbit 

numerical weighting factors, determined by preliminary 
MVS runs 


Minimization of the above cost function provides a fuel optimal trajectory to a 
final orbit which is very nearly synchronous and equatorial. The parameters 
which are allowed to be varied in achieving the above minimization are the thrust 
vector angles a and B (defined in Fig. 3) during the burns and the start times 
of each burn. 


In order to achieve convergence it was necessary to restrict somewhat the scope 
of these variables. The angle between the velocity vector and the component 
of the thrust vector in the orbital plane, a, was assumed to be a linear func- 
tion of time during each burning period. In this way the time history of ot 
during each burn could be represented in the optimization by two parameters, the 
value of a at the beginning of the burn and the value of a at the end of the 
burn. To determine a at any intermediate point the program used linear inter- 
polation. The out of plane angle of the thrust vector, B, was assumed to be a 
constant for each bum. In fact, for low altitude bums the value of p was 

where AV = /"(thrust acceleration) dt 
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V = Velocity Vector 

R = Radius Vector 

N = Out of Plane Vector 

T = Thrust Vector 

a = Angle Between Velocity 
Vector and Projection of 
Thrust Vector into Orbital 
Plane . Positive Outwards . 

0 = Out of Plane Angle of 
Thrust Vector. Negative 
in the Direction of N. 


Figure 3* 


Definition of Thrust Angles 



specified as 7*8 cleg, since this value was found to produce the desired 
2 deg and 26.5 deg plane changes during low altitude and high altitude 
burns respectively. This plane change split is the optimal way to perform 
the required total 28.5 deg plane change with impulsive burns and was found to 
be optimal as well for these low acceleration trajectories by simulating other 
plane splits about this value. The value of p during the final circularization 
burn was iterated upon by the MVS Program, as were the start times of the bums . 
The above restrictions on the variables allowed quick convergence of the itera- 
tive scheme without significantly compromising the solution. 
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3* RESULTS OF THE STUDY 


3 »1 Two Bu rn As cent Mission 

For the two burn ascent to synchronous-equatorial orbit the simulation 
was developed so that cutoff of the first burn was achieved when apogee 
reached synchronous altitude and cutoff of the second burn when the Tug's 
velocity reached synchronous velocity. Start times of the burns and 
thrust angles were iteratively determined as explained earlier. The 
resulting trajectories for both the storable and cryogenic Tugs are shown 
in Tables 2 and 3, respectively. The first trajectory shown in each Table 
is the pseudo optimized trajectory obtained when the in-plane thrust 
during the first bum is assumed to lie along the velocity vector (tangen- 
tial in- plane thrust, = 0 deg). The second trajectory is the optimized 
trajectory obtained when or^ is not constrained to be zero. It is clear 
from comparing total AV figures that the tangential in-plane thrusting 
during the first burn is an excellent approximation to optimal thrusting 
as far as overall AV is concerned. 

The higher gravity losses associated with lower thrust- to- weight ratios 
is evidenced by the fact that the Ml, Wo U (10,000 lbs) thrust cryogenic 
Tug required about 69 mps (225 fps) less to complete the ascent than does 
the 28, 912N (6,500 lbs) thrust storable Tug. Likewise, the burn times for 
the cryogenic Tug are considerably shorter as is the total elapsed ascent 
time. For comparison, the average impulsive thrusting solution is 4,235 mps 
(13>895 fps), so that overall gravity losses are l6l mps (527 fps) for the 
storable Tug and 91 m PS (299 fps) for the cryogenic Tug. 

It is interesting to note that, when the program is allowed to optimize q^, 
the resulting increase in perigee altitude of the transfer orbit is 

Complete optimization was limited by the assumptions on the form of 
the control. 
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less than when Is constrained to be zero. In all cases the inclina- 
tion of the transfer orbit is quite close to the desired 26.5 deg. 
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Table 2. Storable Tug 

Two Burn Ascent to Synchronous- Equatorial Orbit 


Thrust 
I Level 
(lbs) 


Initial 

Orbit 


First 

Burn 


Transfer 

Orbit 


Second 

Bum 


Final 

Orbit 


6,500 


278 X 278 kra 
(150 x 150 n mi) 


i = 28.5 deg 


^ = 0 deg 


6 1 = 7*8- deg 


, v _ 2 , 6 36 mps 

AV - (8,648.8 fps) 

At = 30m 56s 


-2.1+9 

to -,?2 deg 


906 x 35,787 Ion 
(489 x 19,323 n mi) 


i = 26.47 deg 


0 O = -37*1 deg 


AV 


1,763 KPS 
(5,785.2 fps) 


At = 10m 31s 


35,779 x 35,779 tan 
(19,319 x 19,319 n mi) 

i = - .0001 deg 


elapsed time = 

5h 38m 19's 

Total AV = 4,400 mps 
(14,434.0 fps) 


6,500 


278 x 278 km 
(150 x 150 n mi) 


i = 28.5 deg 


<y = -9.21 to 
4.02 deg 


0, = 7-8 deg 


AV - 2 > 62T mps 

AV " (8,618.7 fps) 
At = 30m 51s 




- 2.56 to 
- *73 deg 


828 x 35,787 km 
(447 x 19,323 n mi) 


i = 26.5 deg 


fU = -36.9 deg 


AV 


1,769 mps 


(5,803.6 fps) 
At = 10m 35s 


35,776 x 35,779 km 
(19,317 x 19,319 n mi) 

i = .003 deg 


elapsed time = 

5h 37^ 4ls 

Total AV » 4,396 mps 
(14,422.3 fps) 














Table 3- Cryogenic Tug 

Two Burn Ascent to Synchronous-Equatorial Orbit 


Thrust 

Level 

(lbs) 

Initial 

Orbit 

First 

Burn 

Transfer 

Orbit 

Second 

Burn 

Final 

Orbit 



= 0 deg 


a P = -1.75 

to -.53 deg 

35,857 x 35,859 tan 

(19,361 x 19,362 n mi) 


296 x 296 km 
(l60 x loO n mi) 

6 1 = 7.8 deg 

617 x 35*863 km 
(333 x 19 , 3 ^ n mi) 

?2 = -36.2 deg 

i = .003 deg 

10,000 

i = 28*5 deg 

*v - 2 ?5^8 mps 
AV " ( 8 , 358.0 fps) 

i = 26.42 deg 

. v _ 1,779 mps 
" (5,836.4 fps) 

elapsed time = 
5h 31m 53s 



At = 21m 37 b 


At = 9m 23s 

Total AV = 4,327 mps 
(14,194.4 fps) 



ol = -10.0 to 
-1.78 deg 


= - 2.39 

to - .82 deg 

35,844 x 35,863 tan 
(19,354 x 19,364 n mi) 


296 x 296 km 
(l 60 x 160 n mi) 

= 7.8 deg 

593 X 35,863 tar. 

(320 x 19,364 r. mi) 

B p = -3o.2 deg 

i = .003 aeg 

10,000 

\ 

i ^ 28.5 deg 

AV _ 2 ,545 mps 
A " ( 8 , 350.0 fps) 

i = 26.45 de b 

. v _ 1,782 mps 
A ' (5,345.2 fps) 

elapsed time - 
5h 30m 55s 



At = 21m 3os 


it = 9m 24s 

■ 

Total AV = 4,327 mps 
(14,195.2 fps) 











3.2 Three Burn Ascent Mission 


The three bum ascent to synchronous orbit (shown in Fig. 2) offers the 
advantages of reduced AV, due to lower gravity losses, and the possi- 
bility of using the intermediate orbit as a phasing orbit to allow the 
Tug to inject into synchronous orbit at any longitude. For this analysis 
several intermediate orbit apogee altitudes were chosen for both the 
cryogenic and storable Space Tugs and optimal low thrust trajectories 
utilizing these intermediate orbits were developed. The results are 
shown in Table 4- for the storable Tug and Table 5 for the cryo- 
genic Tug, For these three bum simulations, both and qv, were 
constrained to be zero. This was necessary to keep computer time from 
becoming prohibitive due to the additional variables associated with 
the third burn. As was noted in the previous subsection, the AV penalty 
associated with tangential in-plane burning instead of true optimal 
burning is small for low altitude bums, so that the trajectories of 
Table 4- and 5 are very nearly optimal. 

Plotting AV as a function of intermediate orbit apogee altitude for the 
storable and cryogenic Tugs yields the curves shown in Fig. 4-. 

For the storable Tug the optimal intermediate orbit apogee altitude is 
approximately 8,519 km (4-600 n mi) and the associated minimum AV is 
4-, 302 mps (l4-,115 fps). The optimal intermediate orbit apogee altitude 
for the cryogenic Tug is approximately 6,4-82 km (3,5°0 n mi) and the 
minimum AV is about 4-, 269 m P s (l4-,005 fps). Use of an optimal three bum 
ascent has thus reduced gravity losses to 6T mps (220 fps) and 33-5 mps 
(110 fps) for the storable and cryogenic Tugs, respectively. 

The dotted lines in Fig. 4- cover the range of intermediate orbits about 
the minimum AV point whose orbital periods differ by as much as 90 minutes, 
which is roughly the period of the Shuttle park orbit. 
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Table k. Storable Tug 

Three Burn Ascent to Synchronous-Equatorial Orbit 





Final 

Orbit 


or, = -2.30 to 35,779 * 35,729 km 

* -<59 deg (19,319 x 19,319 n mi) 

R 3 = -36.84 1 = .005 deg 


_ 1,769 nP s elapsed time = 

~ (5,803*9 fps) 7h 07 m 56s 


Third 

Burn 


At '= 10m 37s 


Total AV - 4,380 mps 
(14,370.6 fps) 


-35-87 deg 


35,766 * 35,792 lor. 
(19,312 x 19,326 n ml) 

1 = .006 deg 


1,789 mps elapsed time 


(5,869.3 fps) 


At = 11m 02s 


8b 36m 23s 

Total AV = 4,302 mps 
(14, 114.5 fps) 


-2.51 to 35,770 x 35,783 km 

* .66 deg (19,314 x 19,321 n 

- 36.3 deg i = .0004 deg 


1,780 mps elapsed time = 

(5,839*8 fps) lOh lim 45s 

10m 51s Total fiV = 4,336 mps 

(14,226.6 fps) 















Table 5. Cryogenic Tug 

Three Burn Ascent to Synchronous- Equatorial Orbit 


Thrust 

Level 

(lbs) 


Initial 

Orbit 


First 

Burr; 


Intermediate 

Orbit 


Second 

Bum 


Transfer 


Third 

Burn 


Final 

Orbit 


10,000 


296 x 296 km 

(160 x 160 n mi) 


i = 28*5 deg 


0 ^ = 0 deg 


* 7.8 deg 


_ 1,044 mps , 
AV ' (3>23.7 fps)| 

At - 10m 19s 


332 x 5,556 km 
(17? x 3,000 r. mi) 


i = 27*5 dec 


0 2 = 7.8 deg 


v 1,432 mps 
A (4,699-4 fps) 

At *= 10m 50s 


378 x 35 ,Sr 3 fcr- 
(204 x 1 ?, 36 - a mi) 


i - 26.4 dec 


ar» ** -1.77 to 
- .50 deg 

63 = - 35-5 deg 


av = mps 

6V (5,683.1 fps) 

At = sin 36 s 


35,857 x 35,859 tan 
(19,361 x 19,362 n mi )| 

i = .017 deg 


elapsed time = 

7h 53m 14s 

Total AV = 4,269 mps 
(14,006.2 fps) 


10,000 


296 x 296 km 

(160 x 160 n mi) 


i = 28.5 deg 


* x = 7.8 deg 


_ 1,388. mps , 

A " (4,552.3 fps) 
At = 13m l4s 


372 x 8,521 km 
(201 x 4,601 n mi) 


i = 27.3 deg 


Qfg - 0 deg 


e 2 = 7.8 deg 


AV = 


1,092 mps 
(3,583.3 fps) 


393 x 35,86: km 

(212 x 19, 36 1 - zl) 


i a 26.4 deg 


^ = -1.78 to 
^ -.48 deg 


0 3 = -35.6 deg 


AV 


1,793 mps 

(5,881.2 fps) 


At - 7 m 57 s 


At = 9ir. 35s 


35,857 x 35,859 tan 
(19,361 X 19,362 n mi)| 

i - .0016 deg 


e Lapsed time = 

8h 31m 44s 

Total AV = 4,272 mps 
(14,016.8 fps) 


10,000 


296 x 296 km 
(l 60 x 160 n mi) 


i ~ 28.5 deg 


B x = 7.8 deg 


_ 1,727 mps 
AV - (5,655.3 fps) 


At = 15m 54s 


426 x 12,596 ter- 
(230 x 6,301 n zi) 


i = 27*0 deg 


° 2 = 0 deg 
B £ = 7.8 deg 


v _ 771 mps 
* - (2,529 fps) 


At = 5m 24s 


443 x 35,863 km 
(239 x if, 364 . =1) 


i = 26.4 deg 


= -1.77 to 
* -.49 deg 

° ^ -35-7 deg 


AV 


_ 1,789 mps 


(5,870.8 fps) 
At = 9® 32 s 


35,857 x 35,859 km 
(19,361 x 19,362 n mi) 

i - .008 deg 


elapsed time 
9 h 2> 12s 

Total AV = 4,284 mps 
( 14 , 055.4 fps) 











Total Required AY (fps) 



Figure 4 


Total AV Versus Intermediate Orbit Apogee 
for Three Burn Ascent Mission 







Thus by choosing an intermediate orbit within the span of the dotted lines, 
the phasing of the Tug can be changed anywhere from 0 to $0 min or, equiva- 

i 

lently, anywhere from 0 to 1 revolution in the park orbit. 

Consequently, by waiting in the park orbit for an integral number of 
revolutions (eight revolutions, at most) and then injecting into the 
correct intermediate orbit within the span of the dotted lines in Fig. 4, 
the Space Tug can inject into synchronous- equatorial orbit at any longi- 
tude. From Fig. 4 it is then clear that to allow for worst case longi- 
tude phasing in a three burn ascent to geosynchronous orbit, a AV of 
4,316 mps (l4,l60 fps) is required for the storable Tug and a AV of 
4,278 mps (14,035 fps) is required for the cryogenic Tug. 

Again the increased gravity loss for lower thrust- to-weight ratios is 
evident. The storable Tug requires roughly 38 mps (125 fps) more to 
accomplish the same mission as the cryogenic Tug. 

The AV required to accomplish the three burn ascent, even allowing for 
longitude phasing, is considerably less than the two burn ascent. For 
the storable Tug a savings of 80 mps (262 fps) is realized and for the 
cryogenic Tug a savings of 49 mps (l60 fps) is realized. 
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4. 


CONCLUSIONS 


Near-optimal three dimensional trajectories from a low earth park orbit 
inclined at 28-5 deg to a synchronous- equatorial mission orbit have been 
developed using the Modularized Vehicle Simulation (MVS) Program for both 
the storable and cryogenic Tug baselines. The finite burn times, due to 
low thrust- to-weight ratios, and the associated gravity losses are accounted 
for in the trajectory simulation and optimization. 

A two burn ascent, employing one burn to depart the park orbit and one burn 
to enter geosynchronous orbit was found to require 4,396 mps ( 14,422 fps) 
for the storable Tug and 4,327 mps ( 14,195 fps) for the cryogenic Tug. 

A three burn ascent mission was also investigated. Here the first burn 
produces an intermediate orbit with a specified apogee altitude. The 
second burn, which occurs about perigee, injects the Tug on a transfer 
orbit to synchronous altitude. The final burn circularizes the Tug into 
geosynchronous orbit. For the storable Tug, the optimal intermediate orbit 
apogee was found to be about 8,519 km (4,600 n mi) and the associated mini- 
mum AV was 4,302 mps (14,115 fps). For the cryogenic Tug an intermediate 
orbit apogee altitude of 6,482 km (3,500 n mi) and AV of 4,269 m P s (14,005 fps) 
was optimal. As a comparison, the ideal impulsive thrust solution, in which 
there are no gravity losses, is 4,235 mps (13,895 fps). 

The Tug can inject into geosynchronous orbit at any longitude if the inter- 
mediate orbit is treated as an ascent phasing orbit. In this case a range of 
intermediate orbit apogees must be allowed and the A V is correspondingly higher. 
The storable Tug requires a range of intermediate orbit apogee altitudes from 
5,556 to 12,686 km (3,000 to 6,850 n mi) and a worst case AV of 4,316 mps 
(l4,l60 fps). The cryogenic Tug requires a range of apogee altitudes from 
3,889 to 10,556 km (2,100 to 5,700 n mi) and a worst case AV of 4,278 mps 
(l4,035 fps). Even with these increases in AV to allow longitude phasing, 
the three burn ascent offers a significant AV savings over the two burn 
ascent mission. 
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